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^L( Abstract 

<~i , We investigate how non-standard neutrino interactions (NSIs) with matter can be gen- 

erated by new physics beyond the Standard Model (SM) and analyse the constraints on 

»>f-N [ the NSIs in these SM extensions. We focus on tree-level realisations of lepton number 

^ ' conserving dimension 6 and 8 operators which do not induce new interactions of four 

CO , charged fermions (since these are already quite constrained) and discard the possibility of 

^^ ' cancellations between diagrams with different messenger particles to circumvent experi- 

mental constraints. The cases studied include classes of dimension 8 operators which are 
often referred to as examples for ways to generate large NSIs with matter. We find that, 
in the considered scenarios, the NSIs with matter are considerably more constrained than 
often assumed in phenomcnological studies, at least O(10~^). The constraints on the 

f~^ \ flavour-conserving NSIs turn out to be even stronger than the ones for operators which 

also produce interactions of four charged fermions at the same level. Furthermore, we 
find that in all studied cases the generation of NSIs with matter also gives rise to NSIs 

^ ' at the source and/or detector of a possible future Neutrino Factory. 
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1 Introduction 

With the near start of the LHC, particle physics will enter a new era. With unprecedented 
energy reach and luminosity, the LHC will allow to clarify the origin of electroweak symmetry 
breaking and look for new physics at TeV energies. In addition, complementary to the LHC, 
future precision experiments in the neutrino sector are aiming at measuring the remaining 
unknown parameters in the lepton sector, i.e. the neutrino mass scale, leptonic CP violation 
and the remaining unknown leptonic mixing angle ^13 [1]. 

Precision neutrino oscillation experiments, for example, are also sensitive to new physics 
beyond the Standard Model (SM). This sensitivity is at the same time a chance and a potential 
problem: On the one hand, there is the chance that such experiments discover new physics, for 
example new interactions of neutrinos at the source, detector or with matter, a possible non- 
unitary leptonic mixing matrix, or even a violation of fundamental principles such as CPT 
invariance or locality. On the other hand, new physics may also lead to confusions of effects 
from new (CP violating) interactions with the leptonic Dirac CP phase, in the standard 
parameterisation of the leptonic mixing matrix, or with the small mixing parameter ^13. 
To avoid such confusion when measuring the remaining unknown parameters in the lepton 
sector, a better knowledge of the constraints on the new physics relevant to these experiments 
is highly desirable. 

With respect to their effects on neutrino oscillations, one convenient way to describe 
new interactions with neutrinos in the electroweak (EW) broken phase are the so-called NSI 
parameters for non-standard neutrino interactions at the source (s^a), detector (e^a) [2] and 
with matter (e^g) [3-5]. They give the relative strength of these interactions with respect to 
the Fermi constant Gp. Among these parameters, the NSIs with matter are comparatively 
weakly constrained, i.e. some bounds on them are even 0(1) [4]. In many analyses, large 
non-standard matter effects are therefore included, whereas possible new interactions at the 
source or detector are set to zero. 

In this study, we investigate how non-standard neutrino interactions (NSIs) with matter 
[6] can be induced by new physics beyond the SM. We restrict our analysis to tree-level 
realisations of lepton number conserving dimension 6 and 8 operators which do not induce 
new interactions of four charged fermions (since these are already quite constrained), and 
discard the possibility of cancellations between diagrams with different messenger particles to 
circumvent constraints. The cases studied include classes of dimension 8 operators which are 
often referred to as examples for ways to generate large NSIs with matter without generating 
NSIs at the source and detector of a neutrino oscillation experiment [4,7]. The goal of this 
study is to investigate the constraints on the NSI parameters if these operators are generated 
by explicit new physics beyond the SM. 

The paper is organised as follows: In section 2 we define the NSIs e^^ with matter 
and the related quantities e^^ which affect neutrino oscillations in matter. In section 3 we 
discuss the possible approaches to realising large NSIs with matter as well as the restrictions 
we impose on our analysis. The generation of matter NSIs from dimension 6 operators is 
discussed in section 4, and updated and improved constraints on the NSI parameters are 
derived. In section 5 we investigate NSIs with matter from dimension 8 operators and derive 
the corresponding constraints. Section 6 contains a summary and our conclusions. 



2 NSIs with matter 

Compared to the bounds on NSIs at the source and detector, the NSIs which can modify 
matter effects are often assumed to be only very weakly constrained. In the following we 
will therefore mainly restrict ourselves to this class of NSIs. The (lepton number conserving) 
NSI four-fermion operators of interest are contained in the following Lagrangian after EW 
symmetry breaking, 



C^si = 2\/2Gi. J]e™/4i>La7'i^J [/l7,5/l 



/ 

+2V2GF Y, e^/^ [i^Lal'i^Lp] [/r75/r] • (2.1) 
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The fermions / which the neutrinos couple to are either electrons e, up quarks u or down 
quarks d, and may be left- or right-handed, a, /3 = 1,2,3 are family indices. Constraints on 
the parameters e^g have been derived in [4,8]. 

Neutrino oscillations in the presence of non-standard matter effects can be described by 
an effective square mass matrix which can be parameterised as 

M2j = [/pMNS • diag(0. Ami,, Amii) • ^^pmns + 2^^(diag(l, 0, 0) + e^^) , (2.2) 

where V = \/2Gi?ne, with Ue being the electron number density. The parameters e™^ are 
given by 

^"i,f _ mji^ mjYi =m _ m,e r,m,u m,d f^n , m,u , r, m,d\ /r, r>\ 
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where n„ is the neutron number density. While the individual £^g and £^g are predicted 
in an explicit extension of the SM, only the combined quantity e'^g is relevant for neutrino 
oscillations in matter. 

3 Strategy to realise large NSIs with matter 

Direct experimental constraints on the effective operators of Eq. ()2.ip can mainly be derived 
through neutrino scattering experiments off electrons or nuclei and are, therefore, rather 
weak. The main goal of this work is to investigate whether these mild bounds can be satu- 
rated in extensions of the Standard Model avoiding stronger constraints from other operators 
generated by the same SM extension as a byproduct. Indeed the operator of Eq. ()2.ip is 
not gauge invariant and its generation in any extension of the SM will usually involve also 
the charged lepton partners of the neutrinos in the SU(2)l doublets. The simplest pos- 
sibility to generate NSIs with matter by SU(2)L-invariant operators would be to promote 
the left-chiral leptons (quarks) in the effective operators of Eq. (j2.ip to lepton (quark) dou- 
blets. However, these SU(2)L-invariant operators also generate interactions of four charged 
fermions. Conservatively estimated constraints on the relevant off-diagonal NSIs with matter 
from such operators range from O(10~^) {e^'^' , ejl^**' ) to 0(10"^) {el^if) for the off-diagonal 
elements [9] - [14], while the constraints for the diagonal elements are rather weak (but still 



stronger than the direct bounds of [4, 8])lj In the fohowing, we wih therefore restrict our- 
selves to extensions of the SM (i.e. to the introduction of additional messenger particles and 
interactions) where no interactions of four charged fermions are generated at tree-level. For 
example, this excludes SU(2)L-triplet fermions or scalars as messengers. Their low energy 
effects include interactions of four charged leptons and the resulting NSIs are known to be 
subject to the above constrains (see e.g. [15]). 

We will furthermore not consider here the possibility that extensions of the Standard 
Model with different messenger particles induce these constrained four charged fermion oper- 
ators but conspire to cancel against each other circumventing the bounds. Such cancellations 
between contributions from different messenger particles with different SM quantum num- 
bers are typically associated with a certain amount of fine-tuning without justification by a 
suitable symmetry. We will consequently disregard the extensions in which the extra parti- 
cles and couplings introduced to produce the NSIs of Eq. (j2.ip also lead to diagrams with 
four charged fermions. More generally, we do not consider the possibility that experimental 
constraints are circumvented by cancellations between contributions from different messenger 
particles. 

The lowest dimensional gauge invariant realisations of the effective interactions of Eq. (j2.ip , 
and presumably the least suppressed ones, are provided by dimension 6 operators. Loop-level 
generation of them will always be suppressed by an additional "loop factor" of 0{1/16tt^) 
and cannot lead to very large NSIs. We will therefore restrict ourselves to tree- level genera- 
tion of effective operators. On the other hand, certain dimension 8 operators [4, 7] containing 
four fermions and two SM Higgs fields are often quoted as examples how to generate large 
NSIs (less constrained than NSIs from dimension 6 operators). The reason why dimension 
8 operators appear promising is that the vev of the Higgs SU(2)l doublets can be used to 
"project out" the neutrino fields from the SU(2)l doublets after electroweak symmetry break- 
ing, thus avoiding to generate interactions with charged leptons instead of neutrinos. In our 
analysis, we therefore also include the class of dimension 8 operators (which contains the 
above-mentioned operators) with external fields L,L, f, f,H,H^ , where / can be /l or f^ 
with /l G {Li,Qi} and f^ € {e^,u^,d^} and where H is the SM Higgs doublet. 

In summary, when we attempt to realise large NSIs with matter, we will restrict our 
search to extensions of the SM satisfying the following restrictions: 

• No new interactions of four charged fermions 

• No cancellations between diagrams with different messenger particles 

• Tree-level generation of the NSIs through dimension 6 and 8 operators 

• Electroweak symmetry breaking is realised via the Higgs mechanism 

In the following, we will always impose these restrictions on our analysis. We will start 
with dimension 6 operators and their generation in section [3] and turn to the class of dimension 
8 operators in section [5l 



^Further relaxation of these bounds (up to a factor of 7) is in principle possible [11, 12], however it would 
require specific arrangements in SU(2)L-breaking. 
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Figure 1: Generation of the anti-symmetric 4-lepton operator by the exchange of virtual singly charged scalars 

Si. 

4 Dimension 6 operator for matter NSIs 

There are only two dimension 6 operators and associated SM extensions satisfying the criteria 
defined in the previous section: the anti-symmetric 4-lepton operator, generated from the 
exchange of virtual singly charged scalar fields (c.f. figured]), and the dimension 6 operator 
modifying the neutrino kinetic terms, generated by the exchange of virtual fermionic singlets 
(c.f. figure [5]). The latter operator generates the NSIs in an indirect way, i.e. after canonical 
normalisation of the neutrino kinetic terms. 

4.1 Constraints on the NSIs from the anti-symmetric dimension 6 operator 

In this subsection we will review and update the bounds on the matter NSIs generated 
from the anti-symmetric dimension 6 operator composed of four lepton doublets (see also: 
[10,16,17]) 

^tir = CpT (^« -L.rL.-L^s), (4.4) 

considering its tree-level generation via singly charged scalar fields Si (c.f. figure [T]), i.e. 
new fields beyond the SM in the representation (1,1,-1) of the SM gauge group G321 = 
SU(3)cxSU(2)lxU(1)y. The dot in Eq. (g^]) denotes the SU(2)l invariant product (where 
indices are contracted with e := io"2). In addition to the SM Lagrangian, we therefore consider 
the additional interaction 

^L = -XapLlia2LpSi + H.c. = Ai^5i(£>Li/^ - TpPLUa) + H.c. (4.5) 

as well as a mass msi for the Si. 

Integrating out the heavy scalars Si generates the dimension 6 operator of Eq. (j4.4p at 
tree level. Written in component fields, it has the form 

^tsr = 4E ^^(^'c.PLi^p)i^,PRm = 2Y,^^Cisrpd.){^,i,PLi^p) . (4.6) 

2 Oi i Oi 

For the coefficients c 'T, f^ we can read off 

ct'j.^^ = -Y ^^ipl . (4.7) 
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Using the definition of Eq. ()2.ip . we find that, for normal matter, only the NSIs 
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(4.8) 



are induced. We note that, since the coupling matrix A^^ is anti-symmetric, the indices a 
and (3 in e^!^^ satisfy a^fi ^ eel 



'a/3 



4.1.1 Bounds from rare lepton decays 

One type of constraints in the above extension of the SM comes from rare radiative lepton 
decays la — > Ip")- Neglecting the masses of the light leptons we obtain 



T{la -^ l/Sj) _ _a_ 
r{la -^ Ipt^ai^fs) 487r 






^—^ m\Gp 



with (5 7^ a, /?. Using the present experimental bounds [18] at 90% confidence level (cl) 

Br{ii -^ 67) < 1.2 • 10^" , 
Br{T -^ 67) < 9.4 • 10"^ , 
Br{T -^ /i7) < 1.6 • 10"^ , 



(4.9) 



(4.10) 
(4.11) 
(4.12) 



together with the experimental values Br{T — > Ut^v^j) = 0.1736 ± 0.0006, Br{T — > Vr^^'e) 
0.1784 lb 0.0006 and Br{ii — > u^j,eVe) ~ 100% [19], we obtain the following constraints: 
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-^ rrii Gp 
mlGp 



< 5.0 • 10""^ 



< 1.0 • 10" 



< 4.4 • 10" 



(4.13) 
(4.14) 
(4.15) 



Comparing them with Eq. (14. 8p we see that only r — > /X7 allows to constrain one of the matter 
NSI parameters. At the 90 % cl this constraint is given by 



_m,eL| 

-/XT I 



< 3.0 • 10" 



(4.16) 



This bound turns out to be comparatively weak compared to the bounds that can be obtained 
from the determination of Gp via /i and r decays under the assumption of unitarity of the 
CKM matrix, as we will now discuss. 



^Analogously, the operator in Eq. (|4.4p also induces the NSIs £^a^^ {0;P 7^ m) ^'Hd £^0^^ io.,l3 7^ t), which 
do not play a role for neutrino oscillations in normal matter. The generalisation of the constraints to these 
NSIs is straightforward. 



4.1.2 Bounds from Gp via fi and r decays and assuming CKM unitarity 

The unitarity constraint on the first row of the CKM matrix is experimentally tested to very 
high precision. The extraction of Vnd is performed through superallowed (3 decays, while Vus 
is measured through kaon decayso In both processes Gp, extracted from /x decays, is used 
as an input. Thus, if we assume that the CKM matrix is unitary, the experimental bounds 
provide excellent constraints on new physics contributions to ^ decays. 

The singly charged scalars Si introduced in Eq. (j4.5p can mediate the decay /x — > eUaVp 
with a / e and (3 ^ ^i. For a = /i and /? = e, in particular, the diagram interferes with the 
SM decay amplitude and the suppression of the process will be linear in each A^^ instead of 
quadratic. At this order in A^^, the Fermi constant extracted from the ^ decay would be 
given by 

G, = GHl + E i^^'l ) = Gf{1 + e^^-) . (4.17) 

Using G^ to extract the values of Vud and V^s from (3 decays and kaon decays leads to 

K7 = TTT^ ' (4-18) 



l + e. 



Hfj. 



where V^g^ denotes the experimentally measured V^d and Vus- Using [19] 

V^^P = 0.97418 ± 0.00027 , (4.19) 

V^^P = 0.2255 ± 0.0019 , (4.20) 

and assuming that the unitarity of the CKM matrix is not affected by the new physics leading 
to the NSIs, we find 

Kr\' + \VuT\' = \^;^ = 0.9997 ± 0.0010 . (4.21) 

(1 + e^/l ) 

Analogous to the case of the fi decay, the decay r — > evv' is modified to 

G.^,,p = Gpil + Y^ ^=%7^) = GFil + e^r^) . (4.22) 

The comparison with the /i decay can now be used to obtain bounds on the universality of 
the weak interactions [19,20], which yields 



(j"r — ''evv \ \- -T E' 



m,eL 

TT 



'^fi—*euu V MM 



1.0004 ± 0.0023 . (4.23) 



^The experimental value of Vub is smaller than the precision of the other two matrix elements in the 
unitarity relation and thus negligible for this discussion. 
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1.9 X 10"^ 




8.4 X 10-^ 



Table 1: Bounds on the NSI parameters ej^^ relevant for neutrino oscillations which are generated from the 
anti-symmetric dimension 6 operator given in Eq. (|4.4p . 



4.1.3 Constraints on NSIs with matter 



Using Eqs. (14.211) and (14.230 and additionally the relation lejlr^^l < \/ £]^jf^ ^Tt^^ derived from 
Eq. (|4.8|) . we obtain the following bounds (at 90 % cl): 



"I, El I 

^l^l I 

m,eL| 
r?i,eL| 

^IT I 



< 
< 
< 



8.2 • 10" 
8.4 • 10' 
1.9-10' 



(4.24) 
(4.25) 
(4.26) 



In summary, the anti-symmetric dimension 6 operator of Eq. (j4.4p can only give rise to 
very specific NSIs (with normal matter), namely to e^^^ , ejlr^^ and e^r^^ . The most relevant 
constraints come from the determination of Gp via /i and r decays (under the assumption of 
unitarity of the CKM matrix). Since only the shown matter NSIs involving el are generated, 
the bounds on e^g (which are defined in Eq. (j2.3p and which are the quantities relevant 
for neutrino oscillations in matter) are the same as the bounds on the corresponding s^f^ 
parameters given in Eqs. (|4.24|) - (j4.26p . The bounds on e™^ are summarised in tabled! 

4.1.4 Additionally generated NSIs at the source and their constraints 

In addition to NSIs with matter, the operator of Eq. (|4.6p can also induce non-standard 
neutrino production at a Neutrino Factory source. The coefficients A^o can mediate the 
decay fi — > euu, coupling an incoming /i with an outgoing Da with a / // and the outgoing 
e with an outgoing up with (5 ^ e. Thus, both the neutrino and the anti-neutrino may have 
non standard flavours. Defining 
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T^Lal I^LIS 



[lels{ 



MJ ' 



(4.27) 



we take into account the possibility of both neutrinos having non-standard flavours. We 
notice that for /3 = //, the NSI parameters e^^ ^ reduce to the e^^ usually considered in the 
literature. Eq. (I4.6p gives: 



~'ea,iiP 
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-" V^Gprnl^ ■ 



(4.28) 



From the bounds of Eqs. (I4.13|) - ()4.15p we can derive the following bounds on the e^^ p: 

(4.29) 
(4.30) 
(4.31) 

As in the case of the NSIs with matter, additional bounds can be derived from the 
determination of Gp through fj, and r decays which allows to derive bounds on the individual 
X^i I'^L^/'^'^iP- Yli I'^e^/'^SiP is constrained through Eq. ()4.2ip and Eq. ()4.23l) can constrain 
J2i l-^er/^^SiP- Similarly ^^ jA^^/m^.p can be constrained by [19,20] 
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G - \ ^^ v^™l Gf 



'1 + E 

which results in the bounds 



1.0002 ± 0.0022 , (4.32) 



V2ml Gf 



\et,,J < 8.2-10-4, (4.33) 



'efi,fj,e\ 



£, 



e/i,^r I 



< 1.8 • 10-^ , (4.34) 



\el,^^,\ < 1.9•10-^ (4.35) 

|e^,,^,| < 5.7-10-3. (4.36) 



4.2 Constraints on the dimension 6 operator contributing to neutrino ki- 
netic terms 

The second possibility to generate matter NSIs satisfying the criteria of section [3] is via the 
dimension 6 operator 

^'Frf = -ctf^'ila • H^) i0{H • Lp) (4.37) 

which induces non-canonical neutrino kinetic terms. After diagonalising and normalising the 
neutrino kinetic terms, a non-unitary lepton mixing matrix is produced from this operator. 
The tree level generation of this operator, avoiding a similar contribution to charged leptons 
that would lead to flavour changing neutral currents, requires the introduction of SM-singlet 
fermions (right-handed neutrinos) which couple to the Higgs and lepton doublets via the 
Yukawa couplings (see e.g. [15]) as shown in figure [21 

^It = -YaiiLa ■ H^)K + H.c. . (4.38) 

When singlet fermions (right-handed neutrinos) with Yukawa couplings and a (Majorana) 
mass matrix are introduced, this can in general lead to two effective operators at tree-level: 
It can, on the one hand, produce the dimension 5 neutrino mass operator (Weinberg opera- 
tor) [25] which generates neutrino masses after EW symmetry breaking and violates lepton 
number. On the other hand, this extension of the SM always leads to the dimension 6 op- 
erator of Eq. (j4.37p which contributes to the kinetic energy of the neutrinos and induces 




Nk 




Figure 2: Generation of the dimension 6 operator contributing to neutrino kinetic terms by the exchange of 
virtual fermionic singlets N^. 

non-unitarity of the leptonic mixing matrix. As we will discuss in section I5.H the constraints 
on the diagonal elements of this dimension 6 operator can be used to constrain the NSIs 
induced by the dimension 8 operator. We note that the dimension 5 (Weinberg) operator for 
neutrino masses does not lead to additional constraints, because it can be suppressed by an 
approximate global U(l) "lepton number" symmetry [15,26-31]. The smallness of neutrino 
masses in this case is not explained by large masses of singlet fermions but by the smallness 
of lepton number breaking effects. Sizable non-unitarity of the leptonic mixing matrix could 
arise from comparatively light singlet fermions (with un-suppressed Yukawa couplings), with- 
out being in conflict with the smallness of the neutrino masses. Generally speaking, both 
operators are not directly related since the dimension 5 neutrino mass operator violates lepton 
number while the dimension 6 operator is lepton number conserving. 

From the point of view of neutrino oscillation experiments, having in mind in particular 
a possible future Neutrino Factory, we will regard right-handed neutrinos with masses Mi 
above a few GeV as "heavy" , such that we can effectively integrate them out of the theorycl 
In the following, we review and update the constraints derived in [21] on the product A^A^^^ 
(where A^ is the non- unitary lepton mixing matrix) for Mi larger than the EW scale Aew; 
and extend the constraints to M, larger than a few GeV but below Aew- In the following, 
without loss of generality, we will always work in the basis where the charged lepton Yukawa 
matrix is diagonal. 

4.2.1 The case Mi above Aeav 

In [21] the diagonal elements of NN' were constrained by the combination of universality 
tests and the invisible decay width of the Z. Notice that, without the inclusion of the invisible 
width of the Z, all the constraints derived would consist of ratios of elements of A^A^^ ^nd 
an uncertainty on their overall scale would remain. This can be easily understood from the 
fact that in the Lagrangian the mixing matrix A^ is always multiplied by the weak coupling 
constant g and, since Gp is measured through the fi decay, the comparison of any leptonic 
process will lead to ratios of the elements of NN"^ . Apart from the invisible width of the Z, 
this can be resolved by comparing leptonic and hadronic processes as in section 14.1.21 Indeed 



^Notice that even if Mi < v it will turn out that \vYai/Mi\ <C 1, such that an effective operator expansion 
is still possible. 

10 



the extraction of the Fermi constant from the fi decay with non-unitary leptonic mixing leads 
to 



G^ = GF^{NN%,{Nm)^^ 



Performing the steps as in section 14.1.21 we obtain 



IC'T + IK^P 



(iViVt) 



0.9997 ±0.0010 



(4.39) 



(4.40) 



fj.fl 



To update the bounds of [21], we replace the bound from the invisible decay width of the Z 
by this more tight constraint. 

Furthermore, the off-diagonal elements of NN' are constrained by rare radiative lepton 
decays, la — > Ip'J- With respect to the bounds derived in [21] we also add here the contribution 
of the diagrams mediated by the heavy right-handed neutrinos. This was not considered 
in [21] where a more model independent approach to the source of non-unitarity (based 
on the so-called Minimal Unitarity Violation scheme (MUV) where an extension of the SM 
by only the dimension 5 Weinberg operator and the dimension 6 operator of Eq. ()4.37p is 
considered) was adopted. Notice also that the constraints on the diagonal elements can be 
used to obtain bounds on the off-diagonal ones when the former are stronger, using: 
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(4.41) 



In combination with the additional constraints considered in [21], we obtain the following 
updated bounds at 90% cl: 



l(iViVt) 



a/3 — Sal3 



'^ I d=6,kini 
Y|Ca/3 I 



^4.0 • 10~3 1.2 • lO^'^ 3.2 • 10"^^ 
< I 1.2-10-'' 1.6-10-3 2.1-10-3 
,3.2-10-3 2.1-10-3 5.3-10-3 



(4.42) 



4.2.2 The case Mi below Aew but above a few GeV 

One might think that the constraints on the generated NSIs could be significantly relaxed if 
the singlet fermions are lighter than Aew- For completeness, we will therefore also discuss the 
situation where Mi is below Aew but above a few GeV. From the point of view of neutrino 
oscillation experiments, right-handed neutrinos below Aew but above the typical energies of 
the experiment can still be considered as heavy (and can thus be effectively integrated out of 
the theory inducing non-unitarity of the leptonic mixing matrix) . In general, the constraints 
on non-unitarity of the leptonic mixing matrix from the decays of particles with masses above 
the masses Mj of the right-handed neutrinos are indeed lost, since all the mass eigenstates 
are now available in the decay and unitarity is restored. Thus, the Z and W decays cannot 
be used anymore, however the constraints on the diagonal elements of A^A^^^ derived from fi 



11 



decays, (3 decays and kaon decays together with the universahty constraints from r and vr 
decays stih apply and these can still be translated into bounds on the off-diagonal elements 
using Eq. (I4.4ip . Only the strong constraint on the e// element from /i — > 67 is lost due to 
the restoration of the GIM mechanism. In summary we obtain the following bounds (at 90% 
cl): 

2 /4.0 • 10^3 1.8 • 10~3 3.2 • lO^-'^^ 

\{NN%p - 5^^\ = ^\cif'''n < I 1.8 • 10-3 1.6 • 10-3 2.1 • lO^^ | . (4.43) 
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4.2.3 NSIs with matter induced by the dimension 6 operator which contributes 
to neutrino kinetic terms 

As discussed above, the dimension 6 operator which contributes to neutrino kinetic terms 
leads to non-unitarity of the leptonic mixing matrix, i.e. to {NN^)ai3 7^ (^q,/3. Therefore (c.f. 
[21]), it gives rise to non-standard matter interactions as well as to non-standard interactions 
at the source and detector, which are related to the matter NSIs. In the following, we will 
review the bounds on the matter NSIs in this case. However, we would like to emphasise 
that the related non-standard interactions at the source and detector may also have strong 
(or even stronger) effects on neutrino oscillation experiments^ 

Using the bounds in Eqs. ()4.42p and (14.430 and taking into account that the interactions 
with the W and Z bosons are modified to Nai and {NN'^)^^, respectively, we can compute 
the bounds on the individual NSI parameters in matter induced by the dimension 6 operator 
of Eq. (j4.37p using the relations: 





= 


2 V 2 °^ ^' + 2 -/^ "'') + 


n 

V2" 


- sin 


J 


V^ d=6,kin 


, (4.44) 
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. 2 Q ""^ d=6,km 

-sm Ow—c^p 














(4.45) 




= 


A ^ . 2n \ V^ d=6,kin 














(4.46) 
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lsin^0^p-;'-\ 














(4.47) 


m,di. 


= 


A 1 ■ 2/, A ^^ d=6,kin 














(4.48) 




= 


1.2/1 ^ d=6,kin 














(4.49) 


Using these relations the parameters e^o defined in 


■ Eq. (I23D 


are 


given 


by (see 


e.g. 


[3]) 


~m _ 
^af3 — 


1 

~2 


/ 2 2 \ 
y Cae ^de + y C^^ Oea 1 + 


2 Ue 


V2' 


d=6,km 






(4.50) 



■'The formalism for a full treatment of neutrino oscillations in the presence of such non-unitarity of the 
leptonic mixing matrix can be found in [21]. The NSI parameterisation of new physics in neutrino oscillations 
can also be applied to the case of non-unitarity. Using the NSI approach takes account of the leading order 
effects of the modified interaction with the W and Z bosons, which are induced by the dimension 6 operator 
in Eq. (|4.37p after EW symmetry breaking and canonically normalising the neutrino kinetic terms. 
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NSIs from c'^=B'fe*" 


upper bound (for Mi > few GeV) 


1 rml 
Peel 


2.0 X 10"^ X 1^-21 


Pe/xl 


9.1 X 10-4 X |^-1| 
(for Mi > Aew: 5.9 x lO'^ x |^ - 1|) 


^er 


1.6 X 10^^ X 1^-1 


\'^^J.^J.\ 


8.2 X 10-"^ X ^ 


^^r 1 


1.0 X 10"^ X ^ 


Ptt 1 


2.6 X 10"^ X ^ 



Table 2: Bounds on the NSI parameters e^p relevant for neutrino oscillations which are generated from the 
dimension 6 operator which contributes to the neutrino kinetic terms, given in Eq. (|4.37p . Values nn/rie for 
the crust and the mantle of the earth can be found in table [31 



Compound 


Crust 


Mantle 


SiOs 


60.6 


46.0 


AI2O3 


15.9 


4.2 


FeO 


6.7 


7.5 


CaO 


6.4 


3.2 


MgO 


4.7 


37.8 


NaaO 


3.1 


0.4 


K2O 


1.8 


0.04 


nn/ne 


1.017 


1.019 



Table 3: Values of the % weight amount of the main constituents of the earth's continental crust [22] and 
mantle [23] together with the mean value of ^^^ inferred from that composition. 



which leads to the constraints 
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[ »^e 
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f H" 


- 1 
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\ne 
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1 Hs. 


- 1 


•2 


\ rip 





d=6,kin\ 



d=6,kini 

d=6,kin\ 
C-er 



1 n 



d=d,kini 



I 2 n, ^' 



I 1 rtn 
I 2 ne 



■MM I 

(i=6,fem| 



1 / nn _ X 



=6,fcm|\ 



I 1 rin „ 
I 2 n, "- 



Q!=6,fcm| 



I iBi 

I 2 n, 



Q!=6,fcm| 
C-7-r 



(4.51) 



/ 



•i.\. v^ d=6,kin 

With ^c^^ 



replaced by their upper bounds given in Eqs. ()4.42p and ()4.43p for Mi above 
or below Aew; respectively. Since the ratio ^^ is in general close to 1, this implies that the 



bounds on |e™ | and |e^ 



-mj 



In 



are significantly stronger than the bounds on the individual e^g 
table [3] the values of the % weight amount of the main constituents of the Earth's continental 
crust [22] and mantle [23] together with the mean value of ^^ inferred from that composition 
are given. Notice that the factor ^^ — 1 means an additional suppression of two orders of 
magnitude of the NSI coefficient. 
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(a) Topology 1 



(b) Topology 2 



I 
(c) Topology 3 




Figure 3: Topologies of tree-level Feynman diagrams which can realise the relevant dimension 8 operators. 
The solid external lines in diagrams (a) and (b) correspond to the fields L, L, f, f, H, H^ (where / can be /l 
or /^ with /l G {Li, Q\} and /r G {ej^, u%^, d^}). In diagram (c) the dashed lines indicate the SM Higgs fields 
H and H^ '. 



4.2.4 Additionally generated NSIs at the source and at the detector and con- 
straints 

A non-unitary neutrino mixing matrix N leads to non-standard interactions at the source and 
at the detector of a neutrino oscillation experiment due to the modified coupling to the W . 
In Ref. [24] it has been shown that, parameterising the non-unitary matrix as N = [l + rj)U 
where 77 is a (small) Hermitian matrix and U is unitary, the NSI coefficients at the source and 
detector can be expressed in terms of -q as e^o = e'^a = rja/s- Since NN"^ = (1 -|- rj)^ ~ 1 -|- 2r), 
the bounds on Eqs. (j4.42p and (j4.43p can be translated into bounds on 



'a/3 



^, = |((iViv^W-M = TC''^" 



(4.52) 



5 Dimension 8 operators for matter NSIs 



We now consider the possibility of generating NSIs with matter from dimension 8 oper- 
ators under the restrictions discussed in section O In particular, we analyse operators 
with L,L, f, f,H,H^ as external fields, where / can be /l or f^ with /l G {Li,Qi} and 
/r ^ i^R' '"r' '^r} including their generation at tree-level. 

To start with, operators with these external fields can be generated at tree-level via the 
three topologies of Feynman diagrams shown in figure [3l Scanning over these topologies 
and possibilities for the external fields to be arranged (up to here there are 51 inequivalent 
diagrams for topology 1, 11 for topology 2 and 3 for topology 3) as well as over the virtual 
states allowed to be interchanged, we find only three classes of possibilities which satisfy the 
criteria of section [Hli We will now discuss these three cases as well as the corresponding 
constraints on the NSIs with matter. 



^We note that in some of the discarded cases the interactions of four charged fermions appear after canon- 
ically normalising the kinetic terms of the fields. 
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5.1 Case I: Coupling fermionic singlets to lepton and Higgs doublets 

One generic possibility to generate NSIs with matter via dimension 8 operators while satisfy- 
ing the criteria of section [3] is to couple two pairs of lepton and Higgs doublets to SM singlet 
fermions (right-handed neutrinos) N^, via the Yukawa interactions of Eq. (|4.38|) . Examples 
for diagrams where the external fields / are either e^ or Li are shown in figure IHa) and (b). 
Similar diagrams exist with quarks as external fields. The resulting dimension 8 operators 
have the form [4, 7] 



C 



NSI 



-a/3 



i'f'^tr 



iL^-H^)rf-iH-Lp) 



(5.53) 



These operators are often quoted as examples how to realise very large non-standard matter 
effects. To summarise the classes of diagrams which realise the dimension 8 operators of 
Eq. ()5.53p . we introduce another effective non-renormalisable operator (c.f. figure Hl^c)), 



/:fi = P^MmffNnj 



-'int 



»J 



(5.54) 



In the examples in figure H] this operator is generated by the exchange of a virtual singly 
charged scalar field or by an inert Higgs doublet (which does not get a vacuum expectation 
value). In the following, we will assume that the right-handed neutrinos are heavier than the 
typical scale of a neutrino experiment, such that they can be effectively integrated out of the 
theory and the dimension 8 operators of Eq. (15.530 remain. 

To derive constraints on the NSIs generated by the operators of Eq. (j5.53p . we first match 
the "full theory" with the dimension 8 operators. This leads to the following relation for the 
coefficients c^ : 



d=8,f,I 



E 



Yf 



I3i if) ^jc 



(5.55) 



The corresponding NSI parameters are given by 



'0/3 I 



2 d=8J 



4^/2G/ 



4^/2 



E 






«j 



]a 



Mi Gf Mi 



(5.56) 



We continue by noting that 



Yf3iP^ 



(/) yt 



U 



ja 



^^ Mi Gf Mi 

ij ■> 



< 



2Gf 
2Gf 



\ 



Yf. 






M 



Y, 



aj 



Mi 



d=6,kinii d=6,kini 



'-(3/3 



(5.57) 



The dimension 8 operators 



where p^^' is the modulus of the largest eigenvalue of p\- . 

of Eq. (15.531) thus turn out to be constrained by the bounds on the dimension 6 operator 

contributing to neutrino kinetic energy. 
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m 



s. 



Ni 



(a) 




L'l 



m 



H' 



^ 



K 




(b) 




(c) 

Figure 4: Diagrams (a) and (b) : Generation of dimension 8 operators which induce NSIs with matter after EW 
symmetry breaking in extensions of the SM with fermionic singlets N^ coupUng to lepton and Higgs doublets 
(via topology 1 of figure [3}. Similar diagrams exist with quarks as external fields. Diagram (c): Feynman 
diagram which summarises the diagrams (a), (b) as well as similar diagrams with quarks, introduced in order to 
simplify the discussion of constraints on the generated NSIs. pi^ are effective operators defined in Eq. (|5.54p . 
/ stand for /l or /^ with /l G {I/i,Qi} and /^ £ {e'k,u'B_,d's,}. 
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Now we can use the constraints from the dimension 6 operator contributing to neutrino 
kinetic energy given in Eq. (|4.42p . This leads to the bounds (at 90% cl) 

6.4.10-4 l.l-lO-n (^) 

5.8-10-4 7.3-10-M^, (5.58) 

7.3 • 10-4 1.9 • 10-7 ^^ 

which apphes to both left- and right-handed fermions. 

The bounds of Eq. (|5.58|) depend on the quantity p^^'/Gp which gives the relative strength 
of the effective coupling of Eq. (15.540 with respect to the Fermi constant. Let us now quantify 
how large this factor could be: The extra particles beyond the SM required to generate this 
effective coupling are scalar fields which contain (after EW symmetry breaking) an electrically 
charged component. Such electrically charged scalars would have been produced by pairs via 
photon or Z exchange in the s channel by LEP if their masses were lower than ~ 70 GeV [19]. 
For NSIs with quarks, br and Li should be replaced by their quark counterparts in figure 
HU^a) and (b). The extra scalar fields required in this case would additionally be colored and 
much more tightly constrained. Taking the 70 GeV bound and, for example, couplings of the 
scalar fields to the right-handed neutrinos of 0(1), we find that p)", jGp ^ 10. The bounds 

on the NSI parameters l^^a | are thus at least 0(10-^). 

Furthermore, we would like to emphasise that although the dimension 8 operator of 
Eq. (j5.53p itself produces only NSIs with matter, the additionally generated dimension 6 
operator also gives rise to non-standard interactions with matter (c.f. discussion in section 
I4.2.3p . and additionally it gives rise to non-standard interactions at the source and detector. 

5.2 Case II: Coupling singly charged scalars to lepton doublets 

Another generic possibility to select one neutrino and one charged lepton (and to avoid 
generating in addition couplings between two charged leptons) is to couple a pair of lepton 
doublets to singly charged scalar fields iSj, as in Eq. ()4.5p . The possible diagrams, which 
generate dimension 8 operators of the form 

C^tku = ctfJ''\L^a ■ Lp){L, . m{H^H) , (5.59) 

are shown in figure [SJ For NSIs from these dimension 8 operators, the same constraints as 
the ones for NSIs from the anti-symmetric dimension 6 operator (c.f. section HTTI) still apply, 
since after EW symmetry breaking the dimension 6 operator is recovered. 

5.3 Case III: Mixed case with singly charged scalars and singlet fermions 

In the mixed case (c.f. diagram shown in figure [6]) with one coupling of a pair of lepton 
doublets to Si and one coupling of a lepton and a Higgs doublet to N^^^, constraints on 
the NSIs can be derived analogously to section 15.11 Furthermore, they are of comparable 



magnitude since the constraints on v\Yai/Mi\ < v\/\caa ' *""! (c.f. Eqs. (I4.42|) and ()4.43p ) and 
v\Xl^/mSi\ < 2.9 • 10"^, v\Xl^/mSi\ < 9.2 • 10"^ (from the results of section HTT]) are of the 
same order. The associated dimension 8 operators are of the form 

^Pshn = C^,f'\H^Ll){Lp . H){L, • L^) . (5.60) 
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(a) 




H^ 



Si 



H' 



Sj 



H 




(b) 




S, 



5, 



H / 
/ 



Sk 

A 




\ 



(c) 



Figure 5: Generation of dimension 8 operators which induce NSIs with matter after EW symmetry breaking 
in extensions of the SM with antisymmetric couplings of singly charged scalars Si to two lepton doublets (via 
topologies 1, 2 and 3 of figure [3|. Effectively, the two couplings to the Higgs field only contribute to the 
propagators of the 5*^ after EW symmetry breaking and consequently the constraints on the NSIs in this case 
are the same as for the anti-symmetric dimension 6 operator in section ITTI 



18 




Ni 



w 



H' 



S, 




Figure 6: Generation of dimension 8 operators which induce NSIs with matter after EW symmetry breaking 
in extensions of the SM with fermionic singlets A'^ coupling to one lepton and one Higgs doublet as well as 
singly charged scalars Si coupling to two lepton doublets (via topology 1 of figure [3} . 



Again, we would like to note that the couplings required for generating the dimension 8 
operators also give rise to dimension 6 operators which themselves produce non-standard 
interactions with matter as well as non-standard interactions at the source and detector. 



6 Conclusions 

We have investigated how non-standard neutrino interactions (NSIs) with matter can be 
induced by new physics beyond the Standard Model (SM) and derived the corresponding 
constraints. One motivation for this study was that while NSIs at the source and detector 
are typically assumed to be strongly constrained, in many phenomenological studies very 
large NSIs with matter {s^g parameters 0(1)) are considered, saturating their weak direct 
bounds. To justify such large NSIs with matter (while escaping the stronger constraints that 
would stem from their charged fermion SU(2)l doublet counterparts) it is often referred to 
specific classes of higher-dimensional operators. The goal of this study was to investigate the 
constraints on the NSI parameters if these operators are generated by explicit new physics 
beyond the SM. 

In our analysis, we have focused on the tree-level realisations of dimension 6 and 8 oper- 
ators which do not induce new interactions of four charged fermions (since these are already 
quite constrained). We have furthermore discarded the possibility of cancellations between 
diagrams with different messenger particles to circumvent constraints. The cases studied in- 
clude the classes of effective higher-dimensional operators mentioned above, which are often 
referred to as examples for ways to generate large NSIs with matter. A discussion of the 
restrictions and limitations of our analysis is given in section O 

Regarding dimension 6 operators there are only two possibilities which satisfy the criteria 
defined in section (3) the anti-symmetric 4-lepton operator of Eq. (j4.4p , generated from the 
exchange of virtual singly charged scalar fields 5*, and the dimension 6 operator of Eq. (j4.37p 
modifying the neutrino kinetic terms, generated by the exchange of virtual fermionic singlets 
N^^. The latter operator generates the NSIs in an indirect way, i.e. after canonical normali- 
sation of the neutrino kinetic terms. For both possibilities we have derived improved bounds 
on the NSI parameters (c.f. table [1] and table [2]). The bounds on the quantities |e™^| are at 
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least 0(10-2). 

We have then analysed the possibility to generate NSIs with matter from dimension 
8 operators. Performing a systematic analysis of tree-level generations of operators with 
external fields L,L, f, f,H,H\ where / can be /l or /^ with /l G {Li,Qi} and f^ G 
{e^,Mp^, dp^}, we found that our criteria of section [3] require that either lepton and Higgs 
doublets couple to SM singlet fermions (right-handed neutrinos) and/or that two lepton 
doublets couple to singly charged scalars fields S^. Using the LEP bounds on the charged 
components of the additional scalar fields required for realising the dimension 8 operators and 
allowing for Yukawa couplings 0(1), the NSI parameters \s'^g\ are constrained to be below 
0(10-2). 

In summary, we have found that in the considered setup (c.f. section [3|), NSIs with matter 
are considerably more constrained than assumed in many phenomenological studies, at least 
O(10"2). In some cases the bounds on the NSIs from sources which have previously been 
regarded as quite unconstrained have turned out to be even stronger than for conventional 
dimension 6 operators with structures LLQQ or LLLL (which induce NSIs as well as inter- 
actions of four charged fermions). We have furthermore found that the generation of NSIs 
with matter always gives rise to additional NSIs at the source and/or detector of a possible 
future Neutrino Factory. These NSIs at the source and detector can, for instance, lead to 
"zero distance" neutrino flavour conversion effects which can be efficiently looked for in near 
detectors at future neutrino oscillation facilities. While NSIs with matter with a strength 
below 0(10-2) will be difficult to observe at currently planned or running experiments, they 
might be observed at envisioned Neutrino Factories or /3-Beam facilities [1] and their possible 
impact on precision measurements of the neutrino parameters cannot yet be ignored. In 
order to determine the possible new physics effects in such high precision neutrino oscillation 
experiments, searches at near detectors in neutrino oscillation experiments, improved data 
from EW precision tests and rare lepton decays as well as the results from the LHC will play 
a crucial role. 
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